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Abstract 

0\ 

The Bloch-Nordsieck model for the parton distribution of hadrons 
in impact parameter space, constructed using soft gluon summation, 
is investigated in detail. Its dependence upon the infrared structure 
Oh. of the strong coupling constant a s is discussed, both for finite as well 

(— | \ as singular, but integrable, a s . The formalism is applied to the pre- 

diction of total proton-proton and proton-antiproton cross-sections, 
where screening, due to soft gluon emission from the initial valence 
quarks, becomes evident. 
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1 Introduction 

In this paper we address some phenomenological implications of the infrared 
behaviour of the strong coupling constant a s [|J. In particular, we examine 
some models for the total proton-proton and proton-antiproton cross-sections 
and show the dependence of the rise with energy of the cross-section upon the 
small k t behaviour of a s , through the mechanism of soft gluon summation. In 
a previous paper 0, soft gluon summation techniques have been applied to 
develop a model for the impact parameter distribution of partons in hadronic 
collisions. According to this model, the distribution in impact parameter 



1 



space (6-distribution) is the Fourier transform of the transverse momentum 
distribution of the colliding parton pair, and is obtained by using the Bloch- 
Norsdieck technique for soft gluon summation, developed some time ago to 
describe hadronic transverse momentum distributions 0, £|, |[ || . This model 
for the b— distribution of partons is used in the context of eikonal models for 
total cross-sections, and in particular in the context of the eikonal mini-jet 
models, where the rise with energy is driven by the jet cross-section calculated 
from QCD. In order to make full use of QCD for this particular problem, it 
is necessary that not only the energy dependence be derived from QCD, 
but also the b— dependence, at least for what concerns the hard part of the 
cross-section : it may otherwise be possible to obscure the difficulties of QCD 
inspired models through various parameters which are still present in it. One 
of the difficulties is that the QCD cross-section rises too fast with energy to 
be able to accomodate both the early rise (around ^fs = 10 4- 20 GeV) and 
the high energy behaviour at \/s > 200 4- 300 GeV and beyond. In some 
mini-jet models the too abrupt rise of the mini-jet cross-section is softened 
by modifying the small x-behaviour of the parton densities. Our alternative 
proposal, discussed in detail in this paper, is to regulate the rise of the cross- 
section through soft gluon emission. 

In Sect. 2 we present a brief description of the eikonal mini-jet model. 
In Sects. 34 6 we shall analyze the structure of the Block- Nordsieck model 
for the b— distribution of partons, first recalling the main features of the 
model, and then studying , analytically as well as numerically, its behaviour 
employing various phenomenological models for the k t —>■ behaviour of the 
strong coupling constant a s . In all cases, we shall compare our results with 
those from a model in which the matter distribution of partons is obtained 
from the electromagnetic form factor of the colliding hadrons. In the last two 
sections, Sects. 7 and 8, we shall study the predictions of the Bloch- Nordsieck 
model for total cross-sections and shall compare our results for proton-proton 
and proton-antiproton collisions with other models and present data. It will 
be shown that the model, with a singular but otherwise integrable behaviour 
of a s , is flexible enough to accomodate both the early rise with energy as 
well as present data from the Tevatron. 
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2 Eikonal mini-jet model for total cross-sections 



Ever since the first observation of the rise of proton-proton total cross-section, 
the suggestion was advanced that such rise was due to the increasing (with 
energy) number of hard collisions taking place among the hadron constituents 
[^]. This ansatz was subsequently quantified by the mini-jet model, which 
proposes to calculate the total inelastic cross-section from the jet cross-section 
obtained from QCD H [J. The unitarized version of the mini-jet model is 
represented by the eikonalized minijet model [JTOj, [TTj, [EJ, in which the total 
cross-section is given by 



a tot = 2 J d 2 b [1 - e- n ^' 2 \ (1) 

with 

n{b,s)=A{b)[a soft + cr jet ) (2) 

and A(b) a function which represents the impact parameter distribution of 
partons in the collision. In its most intuitive formulation, the overlap is 
obtained from the Fourier transform of the electromagnetic form factors T\ 
and JF 2 of the colliding hadrons, i.e. 

A FF (b) = J^J d 2 q^T x {q)^{q) (3) 

The model which uses this overlap function, hereafter called the form factor 
(FF) model, although attractive, is of course not parameter free, as it depends 
on the scale parameters characterizing the form factors. 

The two cross-section a so f t and <x, ei are respectively a non-perturbative 
term and a function of energy obtained by integrating the QCD jet cross- 
section from a minimum p t value, ptmin, to the maximum kinematically al- 
lowed. This quantity increases with energy at fixed ptmin, depending upon 
various QCD controlled quantities like the parton densities, in particular, 
and very strongly, upon the small x-behaviour of the gluon densities. In 
fact, the kinematic lower limit in the x-integration for the jet cross-section is 
given by x m i n = 4p^ min / 's, and it can be as low as 1CT 6 at Tevatron energies. 
With such small x-values, the jet cross-section grows much too rapidly as 
s increases and so does the eikonalized cross-section. In order to apply the 
mini-jet model to data, a screening effect is obtained either using the much 
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less dangerous limit y/x m i n or softening the small-x singularity with a cutoff 
parameter. In this way, the above model can reproduce the energy rise, but 
with some further modifications, notably in A(b). In particular, in order 
to obtain reasonable agreement with the data it is also necessary to modify 
the simple form factor model, by allowing for different values of the scale 
parameters for the low and high energy region. 

Our approach is different. We believe that the function A(b) is not a 
constant in energy and for the hard part of the collisions we have proposed a 
model in which soft gluon emission is responsible for the b-distribution of the 
colliding partons. Since the overall soft gluon emission summation is energy 
dependent, we expect such model can modify and complement the mini-jet 
model description of total cross-sections. 



3 Bloch-Nordsieck formalism in impact pa- 
rameter space 

The Bloch-Nordsieck distribution depends upon the energies of the colliding 
quarks and gluons and is thus, although mildly, energy dependent. In this 
section we shall recapitulate the main features of this model, whose general 
structure was derived in ref.Q. As described, our proposed impact parameter 
space distribution for a pair of partons i and j is given by 

e -h(b;M,A) 

Abn = 27r/M&e-MWA) (4) 

where 



h(b , m, A) = HSt f m u + ^-% - Mm (5) 

7T Jo k t M-JM 2 -k 2 t 



with Cy = 4/3 for a quark-antiquark pair. In eqs.([|,|5]) the hadronic scale M 
accounts for the maximum energy allowed to each single soft gluon emitted in 
the collision. This quantity depends upon the energy of the colliding parton 
pair and, through this, upon the energy of the initial colliding hadrons. The 
main point of our model is that soft gluon emission destroys the collinearity 
of the colliding partons. Let us distinguish now between valence partons and 
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gluons or sea quarks. In first approximation, gluons and sea quarks can be 
considered as having the same non-collinearity as the initial valence quarks 
which emit them during the hadronic collision (a different case will be that of 
the photons, which we shall discuss in a different paper). To leading order we 
can now assume that the impact parameter distribution of all type of parton 
pairs is the same as that of the valence quarks. This approximation is in the 
same spirit as the one for which the impact parameter distribution in the 
form factor model is given by the Fourier transform of the electromagnetic 
form factors, i.e. matter distribution follows charge distribution. 

In the calculation of total cross-sections with the eikonalized mini-jet 
model, the distribution @ appears convoluted with parton densities and 
jet cross-sections. In ref.0, we proposed to write the average number of 
collisions at impact parameter b as 

n(b,s) = n soft (b, s)+J2 J ^p- J -p-fiMfjM J dz J dp 2 t A BN (b, gtwas) ^^ 

(6) 

where fi are the quark densities in the colliding hadrons, q ma x is the maximum 
transverse momentum allowed by kinematics to a single gluon emitted by the 
initial qq pair, z = Sj e t/(sxix 2 ), and is the differential cross-section for 
process 

qq -> jet jet + X (7) 

for a given p t of the produced jets with cm. energy J~§jet- The jet pair 
in process (0) is the one produced through any subprocess initiated by the 
valence quark-antipark pair, thus it could be gluon jets, or quark jets. In 
Fig. 1 we show some typical subprocesses which contribute to ([/]). For high 
energy and low p t , most of the jets are produced through scattering of gluons 
emitted by a valence quark pair which continues undetected after emission. 
In principle, an exact calculation of this model for n(b, s) would require to 
know A(b, q max ) f° r eacn subenergy s of the quark- ant iquark pair because for 
process ((7|) 

W(*) = ^(l-^f) (8) 

and then one would need to calculate n(b, s) for each s value, through convolu- 
tion for all parton densities and all subprocesses. This procedure is at present 
unpractical for this problem, since the ^-parameter dependence applies to the 
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Figure 1: Two typical subprocesses contributing to the rise of the total 
proton-proton cross-section 

initial valence pair. What is available, through various parametrizations, is 
parton densities after Q 2 evolution, for all type of partons, whereas the above 
formulation would require to apply corrections and evolution in expressions 
which depend upon the impact variable b. In any case, before recommending 
to embarque in such a time-consuming integration, one can study the prop- 
erties of the proposed model, adopting some approximations, which allow for 
phenomenological calculations. The approximation described in [TJJ is 



where A BN is the function A BN {b, < q max >) evaluated at the value M =< 
Qmax >> obtained by averaging over all parton densities and jet subprocesses. 
In the next section we shall evaluate < q m ax > for different energies of the 
colliding hadrons and for different ptmin values. 




(9) 
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4 The scale dependence : q max 

Using the expression 

_< > 2 ^ ; ^ fi/a{xi) j ^ fMx2) f{dz) 

(10) 

with z min = 4pt min /(sxiX2), one can plot the quantity M as a function of 
•^/s for different values of Ptmin- This is shown in Fig.|, where we have used 
GRV(LO) [13| parton densities for proton proton collisions. 




" 2 3 4 

10 , 10 10 

Vs (GeV) 

Figure 2: The maximum value kinematically allowed for the transverse mo- 
mentum of the single gluon, averaged over densities (GRV-LO parametriza- 
tion) and for different ptmin values, as a function of the cm. energy of the 
colliding protons. 

One sees that, for yfs ~ 50 4- 10 4 GeV, the range of values for M is 
between 0.5 and 4 GeV for ptmin = 14-2 GeV. For these typical values, one 
can now calculate h(b; M, A) and subsequently A(b, M). Our point of interest 
in this paper is also to relate the rate of rise of the total cross-section with the 
behaviour of a s in the infrared region. The stronger the singularity as k t — > 0, 
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the larger h(b;M,A), the faster Abn goes to zero and the stronger will be 
the suppression produced by soft gluon emission. We shall now quantify this 
statement with numerical calculations. 



5 a s dependence in the function h(b; M, A) 

We start by showing how the b-dependence of h(b; M, A) varies according 
to the behaviour of a s in the very low k t region. Because of the many 
uncertainties we shall work with the one-loop expression for a s and shall use 
two different models, each of them characterized by a set of parameters, i.e. 
the frozen a s model used in [Tj| where 

/»2\ _ 127T 1 

a* ) ~ 33 _ 2Nf ln[(fc 2 + a 2 A 2) /A 2)] ^ 

which depends upon the parameter set {A, a} and in which a s goes to a 
constant value as kt goes to zero. An altogether different model is the singular 
a s model, described in with 

/,9s 127T V 

= m 9AM, n 1 (12) 

which coincides with the usual one-loop expression for large values of k±, 
while going to a singular limit for small k±. In this model, a s depends upon 
the parameter set {A,p}. The singular expression of eq. flT2"|) is inspired by the 
Richardson potential fT6j used in quarkonium spectroscopy. The Richardson 
potential can be connected to a singular l/kj_ behaviour of a s in the infrared 
limit, a singularity which is not dangerous in bound state problems, where the 
Schroedinger equation selects only those solutions for which the momentum 
is fixed by the stability condition. For this problem, and as discussed in 
[0], the expression we have chosen should be considered as a toy model, 
in which the singular behaviour of a s (if any) can be modulated through 
the singularity parameter p. One should also notice that the singular limit 
of the above equation is not an observable. Phenomenologically, one never 
measures a s in the k± — > limit, since this limit corresponds to emission 
of a very soft gluon, in which case summation, and hence integration over 
k±, is mandatory. In other words, what really matters is the integrability 
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of the function, since observable quantities (soft gluons are observed only 
as overall energy momentum imbalance carried away by soft particles, but 
not measured individually) always involve an integration over the infrared 
region. In what follows we shall always use the set a = 2, A = 0.2 GeV for 
the frozen a s model, whereas for the singular case, while we shall vary the 
singularity parameter p, we shall adopt the value A = 0.1 (jeFljlTfl. 

Let us now examine the function h(b; M, A). This function does not allow 
for a closed form expression, and needs to be numerically evaluated. Useful 
analytical approximations can be found in the appendix. 

The dependence of the function h(b; M, A) upon the infrared behaviour 
of a s is shown in Fig|| where we have plotted in the same graph the exactly 
integrated expression for the function h(b, M, A) for the frozen case, eq.(ll), 
and for the singular case, eq.(12), for two values of the parameter p. For each 



50 
45 
40 
35 
.30 
•25 
20 
15 
10 
5 




1 af 9 with p=5/6 (dashes) 



t" < ng with p=3/4 (full) 



a' roz (dots) 




8 10 12 14 16 18 20 



b(GeV" 1 ) 



Figure 3: Comparison between numerically integrated expressions for 
h(b, M, A) for singular and frozen a s . M = 1 (lower) and 4 GeV (upper). 



case, we have evaluated h(b) for the two values M = 1 and 4 GeV, which 
correspond to the interesting range t/s = 50-i-10 4 GeV, for j> t min = 1-^-2 GeV 
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(see Fig. |2[). Although at very small b (b < 0.2 GeV^ 1 ) the values are not 
very different, at larger b- values there is an increasing discrepancy between 
the two formulations. The large b region below m 10 GeV~ x is the one which 
matters most for the total cross-section analysis, where the figures then show 
that the infrared behaviour of a s plays an important role in the rise of the 
cross-section. In the next sections we shall study the difference in A(b) and 
then in the number of collisions, given the same jet-cross-section. 



6 The overlap function A(b) 

In this section, we shall calculate numerically e -fc ( fe;M ' A ) and the normalized 
A(b), for the two cases, frozen and singular a s . We show in Figs.|4],[| the 
normalized function A (b) for the frozen and singular a s possibilities, using, 
for the latter case, three different values of the parameter p which regu- 
lates the singularity. In both figures we also show the comparison with the 
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Figure 4: The overlap function 
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Figure 5: As in Fig.^ for singular 
a s , for various values of the pa- 
rameter p. 



function A(b) in the form factor model, according to which matter density 
in the proton is given by the electromagnetic form factor. With the usual 
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parametrization 




2 



v 2 = 0.71GeV 2 



(13) 



the overlap function A{b) in the form factor model has the expression 



In each figure, the various curves correspond to varying the scale M as 
described in the previous section, so that they include a range of energies 
\fs = 50 4- 10 4 GeV for a range of ptmin between 1 and 2 GeV. We see that 
the frozen a s case is more similar to the form factor model, especially at low 
medium energies, (50 — > 100 GeV), when the proton is not yet exhibiting the 
full QCD behaviour. This is different from the singular case, where the func- 
tion A (b) is always falling with energy more than in the form factor model. 
The more singular a s as k t — > (larger p values), the more concentrated at 
small impact parameter is the overlap function and hence the less important 
the large b-values. This will have as physical consequence that as the cm. 
energy increases, the non-collinearity of the initial state due to soft gluon 
emissions will accordingly increase. Clearly this will signify a much more 
noticeable effect of soft gluon straggling on the total cross-section. We shall 
now see this effect on the average number of collisions n(b, s). 



In the eikonalized mini-jet model, the quantity which contains the energy 
dependence of the total cross-section, is the average number of collisions 
n(b,s). At low cm. energy of the colliding particles, this number is dom- 
inated by contribution from soft, non-perturbative type events, while the 
QCD component, mini-jet like, slowly rises, reaching a comparable size in 
the 200 4 300 GeV region. As mentioned in the first section, one can ap- 
proximate the average number of collisions in the entire region as 




(14) 



7 Average number of collisions 




(15) 



with 




(16) 
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and 

n har d(b, s) = A BN (b; M, A)cr jet (s,p tm in) (17) 

To study the b-behaviour, we shall introduce the soft term, by using the form 
factor model for A(b) described in the previous section, which is consistent 
with a low energy model of the proton in which only valence quarks play a 
role in the scattering. In this model all the energy dependence comes from 
the cross-section term : we will parametrize u so ft so as to reproduce, through 
the eikonal, the low energy behaviour of the total proton-proton and proton 
anti-proton cross-sections. We found, as best fit to the low energy data with 
an eikonal formulation with rihard = 

<f t = ^ + ^ (18) 

and 

< t = 47 + ^T + ^ (19) 

where E is the proton energy in the Laboratory system in GeV and the 
cross-sections are in mb. For <jj et we use GRV(LO) densities to evaluate 
the proton-proton jet cross-section and two different values of Ptmin = 1.2 
and 2 GeV, the latter being the one for which the total cross-section in the 
FF model passes through the CERN data points at y/s = 546 GeV. In a 
subsequent section, when we shall try to fit the total cross-section data, we 
shall use other ptmin values. We can now plot the entire n(b, s) as a function 
of b, for various values of the center of mass energy y/s, which corresponds to 
various values of the scale M, as described in the first section. We show this 
behaviour for the frozen and singular a s case in Figs.||, for p tm ,in = 2 GeV . 

For the frozen a s model, shown in Fig.[], the results are compared to a 
straightforward application of the form factor model, i.e. with 

n FF (b, s) = A FF (b) [a ^ + a jet ] (20) 

We see that at y/s = 100 GeV, there is still no difference betwen the two 
models. On the other hand, as the energy increases, the BN model shows 
a stronger suppression of the large 6-contribution. For the singular case, in 
order to show variations with the singularity parameter p, we plot in Fig.|7| 
the result for different p-values. Notice that the p-dependence is related to 
the values of k t probed, i.e. by the M values, which are smaller the smaller 
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Figure 6: The average number of 
collisions for the frozen and singu- 
lar a s case in comparison with the 
FF model at various cm. energy 
values. 




Figure 7: The average number of 
collisions for the singular a s case 
for different values of the singular- 
ity parameter p. 



Ptmin is. Thus, for ptmin = 2 GeV for instance, there is very little difference 
among the various curves, at any given energy. This reflects the fact that the 
upper integration limit in eq.([5|) is a relatively large M (3 -S- 4 GeV) value, 
so that the overall function is not very sensitive to the infrared region. It 
should be noted that for smaller ptmin values, like the ones actually used for 
fitting the total cross-sections in the next section, the dependence upon p is 
much more noticeable. We show one such case in Fig.|7]. 

The next figure, Fig.[8], shows a break down of the average number of 
collisions into the soft and the hard component. In the present analysis 
we are not changing the soft component, which appears as the dash-dotted 
curve, and the figure shows, at a given high (LHC) energy, how the hard 
part would be different in the three models, i.e. in general more peaked 
at small b for the Bloch-Nordsieck model, and in particular falling faster 
the stronger the singularity of a s . At lower energies, where the mini-jet 
contribution is less important, these discrepancies would be much reduced. 
So, in this picture, while keeping a similar b-distribution at low energy, we 
quantitatively enhance small b-collisions at high energy, though QCD soft 
gluon emission. The change in the 6-distribution introduced in the hard 
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Figure 8: Soft and hard compo- 
nent of n(b,s) in the three models 
described in the text. 




Figure 9: The average number of 
collisions in the form factor model 
and the Bloch Nordsieck model, at 
LHC energy. 



component by the different models for A (b) is responsible for the changed 
shape of n(b, s) between the form factor and the other two models. The direct 
comparison among the three models is shown in Fig.(^) where the average 
number of collisions at y/s = 14 TeV is plotted for a choice of the various 
parameters as indicated. Apart from the change in shape, it can be noticed 
that the frozen a s case corresponds to a behaviour intermediate between the 
form factor model, and the singular a s case. We also see from this figure 
that the range of values of the h— parameter most important to the total 
cross-section calculation changes in the different models. 



8 Total cross-sections 

Before attempting the last relevant phenomenological exercise for the cal- 
culation of the total proton-proton and proton-antiproton cross-section, we 
shall first show how the integrand in eq.(|l|) changes with energy and which 
values of b are most relevant for the calculation of the total cross-section in 
the various models for A{b) we have just described. We must stress that 
this is not an optimization of the many parameters from which this model 
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depends : rather an exercise to show how the Bloch-Nordsieck model for 
the impact parameter distribution affects the total cross-section behaviour 
in the eikonalized minijet model and how the behaviour of a s in the infrared 



region is related to the rise of the total cross-section. This is done in Fig. 10 



The figure shows how much the integrand of eq.([|) is peaked at different 
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Figure 10: The integrand of the eikonal formulation for a tot , for Ptmin = 
2 GeV in the three different models described in the text, for a range of cm. 
energy values 100,1000,10000 GeV 



6-values as the energy increases, but also as the model for A(b) changes. And 
it indicates that the rise with energy of the area under the curve, i.e. the 
cross-section, at the same energy shrinks for the more singular a s behaviour. 
Finally in Fig.ITT] we show the comparison of this model with proton- 
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proton and proton-antiproton data. For proton-proton, we only show data 
up to ISR energies, since the existing data points in the TeV range are 
extrapolations from cosmic ray data JTB| from p — air collisions and are partly 
model dependent [pj| . For the proton-antiproton data, we have plotted all 



the data points published so far from the CERN SppS |2(J and FNAL pl| 
experiments. This introduces a larger band of uncertainty that it is usually 
shown, but the purpose of this paper is to indicate the potentiality of the 
Bloch-Nordsieck model rather than to do a best parameter fit, and we have 
opted for a comparison of our results with the full experimental picture. 

We have studied three different formulations of the eikonal mini-jet model, 
one for the form factor model and two for the soft gluon summation model. 
To choose the parameters of the mini-jet description, we have selected those 
Ptmin values which would ensure that the curve can reach the high energy 
points : for the form factor model this can be accomplished with ptmin — 
2 GeV but, as often stressed, with such value it is not possible to fit the 
early rise of the data. A lower value of ptmin would on the other hand give 
curves which rise too much at higher energy and miss the points. Going to 
the soft gluon summation model, it must be noticed that since this model 
has an energy dependence in the b-behaviour in addition to the one in the 
jet-cross-section (common to all the models), one can expect that a smaller 
Ptmin could be used, thus allowing for the earlier rise. In fact, the high 
energy data, for the frozen ot s model, can be met with ptmin = 1-6 GeV. 
Although with this value, the cross-section starts rising sooner than in the 
form factor model, still it is impossible to fit both the early rise as well 
as the high energy points. This model depends not only upon p t min value, 
but also on the scale a which regulates the infrared behaviour of a s : the 
smaller a, the more singular the behaviour and the easier to fit the early 
rise. Finally, we show the results for the singular a s case, with a particular 
choice of the parameter p which regulates the singularity of a s . We can 
choose now a rather small value of ptmin to reproduce the early rise, since at 
higher energy the increased soft gluon emission reduces drastically the large- 
fa contribution to the cross-section and does not let it rise as much as in the 
other models. Our results are compared with a multiparameter fit from a 
QCD inspired model [0, which has recently been used to successfully predict 
photon photon total cross-sections p2| . These results are not very different 



from the ones obtained using the Regge-Pomeron exchange picture | 23| , but 



the model in p2] is closer in spirit to the one discussed here, with the energy 
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rise due to the rise of the QCD jet cross-section. 

The results of this figure shows that it is possible to have a rise in agree- 
ment both with the intermediate energy data as well as with the Tevatron 
data : this result is obtained using a single eikonal function, usual QCD par- 
ton densities and minijet cross-sections with ptmin i n the 1 -j- 2 GeV range. 
To follow the beginning of the rise, one needs a rather low Ptmin- 111 general 
such low values imply too fast a growth of the total cross-section, in our case 
this fast growth is tampered by the increasing number of soft gluon emission 
phenomena at small k t . 

9 Conclusions 

We have presented a detailed numerical analysis of a Bloch-Nordsieck ap- 
proach to the impact parameter distribution of partons in the context of the 
eikonal mini-jet model for total hadronic cross-sections. We have shown that 
the proposed soft gluon summation expression plays an important role in 
softening the rise of the cross-section due to mini-jets and have studied the 
role which the infrared behaviour of a s plays in it. 
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A Approximate Expressions for h(b, M, A) 

In this section, we show some analytic approximations to the function h(b, M, A). 
Here we shall restrict our attention to values of M relevant to the total cross- 
section calculations, i.e. values in the few GeV range. Since as the total 
cm. energy increases, M increases from 0.5 to 4 GeV, the region of b > 1/M 
(b < 1/M) corresponds to values of b larger (smaller) than 2.5 GeV -1 , at low 
v^s, down to 0.2 GeV^ 1 for the highest y/s values. In other words, in the 
integration, small and large b-values are an energy dependent concept : at 
very small y/s, small b, i.e. b < 1/M means values of b less than 2.5 GeV^ 1 , 
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whereas at very high energy large b-values mean b > 0.2 GeV. We shall now 
start studying h(b] MA) in the frozen a s case, and distinguish three cases : 

1. bM < 1 

2. bM > 1, baA < 1 

3. bM > 1, baA > 1 

In order to obtain a closed form expression to better study the function, 
we shall adopt the following approximations : 



a s (k < aA) = a s 



12tt 



27ln(a 2 ) 



a s (k > aA) 



12tt 



27/n-g 



In 



In 



M-VM 2 -k 2 
M— y/M 2 -fc 2 



2 In 
2 In 



2M 

T" 

M 

T 



k Pi and 

for values not in the infrared region. 



x 



1 — Jo(x) — — X < 1 

1 - J (x) = 1 X > 1 

Then, one can break the integral from — > M into various intervals in which 
one can approximate the integrand and perform the integration. According to 
the three cases indicated above, one then obtains the following approximate 
expression : 



bM < 1 

2c F 



(A.l) 



h(b,M,A) = 



7T 



b 2 r aA , „ , 2M 
( i , s — / kdk In — — 



+ b 



2c F _ 
< a 

7T \ 

b 2 M 2 



b 2 A 2 a 2 r 01 2M 

, l + 21n — 

8 [ aA 

fa 2 A 2 n A 2 
s- - 1 + 2 — K In 

M 2 M 2 





/■M 










+ 




M 


Zi ^ 




T 



kdk—^- 



il^ 2 
A 2 



- U{a 2 
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For bM > 1, one distinguishes between two cases : b larger or smaller than 
1/aA, so that the integral can now be divided as follows : 



for 
h(b, M, A) 



M 



aA 

b 2 r aA 2M -b 2 r- b lnf - 
q,- / kdkln — ho— / kdk — jr + b 

k 4 JaA 



2c f 



7T 

2c F 

7T 

-b 2 A 2 



2 . 

b 2 A 2 a 2 r 



(A.2) 

r dk lnf" 
fc ln-l 



a 2 - 



1 + 2 In 



2M 
aA 

1 «, M 

7577 + 2 ln 7T 
o 2 A 2 A 

' M \m — 
In — In - ln(Mb) 



+ 



o 2 A 2 



- Z*(a 2 ) 



+ 



or as 



for 
h(b, M, A) 



M 



< — <6 
aA 



2c» 



7T 

2c F 

7T 



_ b 2 /■*,„. 2M o _ r«*dk, M T f 
a s — / kdkln— h 2a s / — In — + 6 / 

2 Jo /c J\ k k Ja 



(A.3) 
^ rffclnf 

aA fc In | 



^[l + 21n(2M6)] + 

o 



2a s |ln(M6) ln(aA6) - X - In 2 (aAo)} 



+ 



A In a 



i M 

aA 



The last decomposition is the one to use to study the large b limit, whereas 
the first one corresponds to the small b limit. 

For the singular a s case we adopt similar approximations, except that 
now 



a s (k < N P A) = b (-) 
a s (k > N p A) 



~k 
b 

Ml) 
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l/2p 



where N„ 



is a number of order unity. For p=l/2, indeed N p = 2 



and the two regions, small and large k, coincide with those in the frozen a s 
case with a=2. The expressions for h(b; M, A) in this case become : 



b < 



h{b; M, A) 



1 

M 
2c F 

7T 

2c F 



(A.4) 



b 2 r 

4< A > 2p i 



dk , 2M -6 2 
In — + b— 



k 2p- 



k 



2U-P 



7T [ 8(1 — p) 



W 2 In 



2M 



AiV, 



+ 



-6 2 M 2 



iV 2 A 2 A 2 , 



v 1 ~P 
M 



4 Jan, 



+ 



and for the bM > 1 case one will have the two possibilities, 



for 



1 



M 

h(b, M, A) 



< b< 
2c F 

7T 

2c F 

7T 



1 



N p A 



2 



10 

b 2 A 2 (A" 2 ) 1_p / _ 2 U 



P 



k 2p-l 

2 In 

M 



rfA; , 2M -b 2 r 
b—A 2p / — — In — + b— 

Jo k 2 ^- 1 k 4 Jn v a 



* hdk H + ij" dk,a * 



(A.5) 

a; 



k In 



1 



iV p A 



1 — p 



+ 



N 2 b 2 A 2 -l + 2\n—b 2 A 2 
p A 

' M ln M 
In l n A _ ln ( Mb \ 

. A 



^ W } ~ u{Np) 



+ 



and the other case 



b > 



1 



N p A 



h(b, M, A) 
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> M 
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2c F 

7T 



■ b 2 A 2p 
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2p 



(b 2 A 2 y 



2 ln(M6) 



+ 



2pN t 



2p 



-2 In 



M 
AN 



+ 



V 



+ 



b In 
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1 



\ lnf lnAL 
In ^A- - 1 + P 



In AC, 



This approximation is reasonably accurate, as one can see from Figs. 12 -13 



where we have plotted both the approximate and the exact expressions from 
the above equations for the two different models for a s . 
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Figure 11: Total p-p and pp cross-sections and comparison with various mod- 
els 
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Figure 12: Comparison between 
the approximate and the ac- 
tual numerical integration for 
h(b, M, A) for various values for 
M, in the frozen a s model 




b (GeV ') 



Figure 13: Comparison between 
the approximate and the actual 
numerically computed expression 
for hip, M, A) for the singular a s 
model 
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